Fluctuation conductivity of thin films and nanowires near a parallel-field-tuned 
superconducting quantum phase transition 
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We calculate the fluctuation correction to the normal state conductivity in the vicinity of a quan- 
tum phase transition from a superconducting to normal state, induced by applying a magnetic field 
parallel to a dirty thin film or a nanowire with thickness smaller than the superconducting coherence 
length. We find that at zero temperature, where the correction comes purely from quantum fluctu- 
ations, the positive "Aslamazov-Larkin" contribution, the negative "density of states" contribution, 
and the "Maki-Thompson" interference contribution, are all of the same order and the total cor- 
rection is negative. Further we show that based on how the quantum critical point is approached, 
there are three regimes that show different temperature and field dependencies which should be 
experimentally accessible. 

PACS numbers: 



The high interest in the physics of a quantum critical 
point (QCP) is motivated by the explosive growth of its 
experimental realizations, driving, in turn, a quest for 
theoretical understanding of quantum phase transitions 
(QPT)Hl- The challenge for theory is two fold: First, it 
is desirable to identify experimentally accessible systems 
experiencing QPT. Second, these systems should allow 
for a systematic and comprehensive theoretical descrip- 
tion. In this respect, the QCP realized in dirty super- 
conductors of reduced dimensions under application 
of magnetic field represents an exemplarily controllable 
system allowing complete experimental exploration of the 
vicinity of the critical point on the one hand and system- 
atic theoretical study on the other. 

In this letter we present a full systematic investigation 
of fluctuation corrections^, to the normal state con- 
ductivity of a thin wire or a thin film in the vicinity of 
a QPT from superconducting to normal state, induced 
by an applied magnetic field. We consider a thin wire 
(or thin film) of diameter d (or thickness t for the film) 
much smaller than the superconducting coherence length 
£. The magnetic field H is applied along the wire (or par- 
allel to the film) and can be parameterized completely by 
a scalar pair-breaking parameter 0] (similar to the effect 
of paramagnetic impuritiesQ) 
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where D is the diffusion coefficient in the bulk sample 8] . 
The superconducting critical temperature T c is related to 
a (see Fig. 1) via the standard equation 



ln(T c /T c0 ) = ^(1/2) - V(l/2 + «/2ttT c ), 



(2) 



where tp is the digamma function and T c q = T C (H = 0) 
is the critical temperature in the absence of a magnetic 
field. One can immediately see that at T — the super- 
conductivity is destroyed at a c o = a c (T = 0) = ttT c o/2j, 
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FIG. 1: Phase diagram of a superconducting nanowire (thin 
film) in a parallel magnetic field (parametrized by a). The 
boundary between the quantum and the intermediate regime, 
shown by a dotted line is different for a wire and film as 
marked in the figure while all other boundaries coincide. 



where I117 k 0.577 is the Euler constant. Realization 
of the superconducting QCP in a superconductor with 
paramagnetic impurities was first suggested in the work 
of Ref. jsl (see also Ref.[lfl|). Since in our case the depair- 
ing parameter a depends on H , it allows for a well con- 
trolled exploration of the QPT and its vicinity by varying 
the applied magnetic field. 

Our main result for the fluctuation correction to the 
normal state conductivity can be presented as 



Sa(a, T) = Sao(a) + Sa T (a, T), 



(3) 



where 8<jq{ol) and 5aT{u,T) are the zero temperature 
and finite temperature contributions, respectively. The 
zero temperature correction is given by 
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where a q = a - 



ird(2 + d) J (2Tr) d a\ \n{a q /a c ) ' 
Dq 2 1 2 and d =1,2 corresponds to a wire 
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and film respectively. Its magnitude decreases monoton- 
ically with increasing field; this leads to a negative mag- 
netoresistance. Note that a negative magnetoresistancc 
was found also in granular superconductors and in 
thin films in perpendicular magnetic field |l2j. Shown in 
Fig. 2 are plots of the dimensionless correction 
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When expanded around the QCP, the dimensionless con- 
ductivity close to the QCP is given by 



Sa (a) = 6a (a = a c0 ) + b(a - a c0 )/a c0 , 



(6) 



with the numerical coefficient b = 0.386 and b = 0.070 for 
a wire and film, respectively. Thus the zero temperature 
critical exponent is 1 for both d = 1,2, at least within 
the perturbation theory. 

In the vicinity of the QCP (T, a — a c <Ca c o ), the field 
dependence of SaT(o:,T) turns out to be more singular 
than that of Sao(a) and for T > a — a c its leading term 
is given by 



wire 
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while for T < a — a r it is 
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The contributions Sa T and Sao become comparable at 



T = T (a) 
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The key point is that the behavior of the fluctuation cor- 
rections to the conductivity depends on the way one ap- 
proaches the QCP and we can identify three regimes in 
the vicinity of the QCP that show qualitatively different 
behaviors as illustrated in Fig. 1. There is a "classical" 
regime for T > a — a c , where the correction is given by 
Eq.0; an "intermediate" regime for a — a c > T > To (a), 
where the correction behaves according to Eq. JHj ) ; and 
a "quantum" regime for To (a) > T where the behav- 
ior crosses over to an essentially zero-temperature-like 
behavior which is not singular and almost temperature 
independent with the fluctuation correction dominated 
by Sao as given by Eas. (|4l5l6f> . 

Since in the quantum region the correction to conduc- 
tivity is negative whereas in the classical and interme- 
diate region it is positive, we predict a non-monotonic 
behavior of the resistivity as a function of the magnetic 
field at finite temperature. The corresponding plots are 
shown in the insert of Fig. 2 for nanowires and thin 
films. Such behavior was indeed observed in experiments 
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FIG. 2: Dimensionless zero temperature fluctuation con- 
ductivity correction© as a function of depairing parameter 
a. Insert shows dependence of the resistivity on the mag- 
netic filed for temperatures T/q c o = 0.01 (upper curves) and 
T/a c o =0.1 (lower curves). Resistivity is normalized to the 
high temperature resistivity R® while the magnetic field is 
normalized to the critical filed H c o at T = 0. 



on amorphous thin films |2j, while for nanowires, to the 
best of our knowledge, it was not reported yet. 

Our approach is based on the diagrammatic pertur- 
bation theory. Note that the technique using the time- 
dependent Ginzburg-Landau formalism that was adapted 
in Refs. 0, , accounts only for the direct "Aslamazov- 
Larkin" (AL) type of contribution^^] to the fluctuation 
conductivity that comes from the charge transfer via fluc- 
tuating Cooper pairs, but misses the zero-temperature 
contribution to the correction. On the contrary, our 
approach, takes care of all the contributions including 
the "density of states" (DOS) part resulting from the re- 
duction of the normal single-electron density of states at 
the Fermi level, and the more indirect "Maki-Thompson" 
(MT) interference contribution [l3. fl5| . At zero temper- 
ature, where the correction comes purely from quantum 
fluctuations, these turn out to be of the same order as 
the AL contribution. 

To derive our main results, we carry out a microscopic 
calculation within the standard framework of temper- 
ature diag rammatic technique in a disordered electron 
system [lfillTlllsj ] in the diffusive limit (inverse mean free 
time r^ 1 ^> T,a). Since we choose d,t < £ and param- 
eterize the effect of field by a as defined by Eq. ( , the 
problem effectively becomes one and two dimensional for 
a wire and film, respectively. The main building block of 
the diagrammatic technique in the presence of the BCS 
interaction is the so-called "Cooperon" , the ladder dia- 
gram that describes coherent scattering by impurities in 
the particle-particle channel. In the presence of a parallel 
magnetic field it is given by 



C(wi,u}2,q) 



2llv6( — LOlUJ2) 

|wi - lo 2 \ + Dq 2 + 2a 



(10) 



where 9 is the step function, v is the density of states at 
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FIG. 3: Diagrams for the fluctuation conductivity: a) 
Aslamazov-Larkin diagram, b)-c) Maki- Thompson diagrams, 
d)-e) density of states diagrams. Full line stands for the dis- 
order averaged normal state Green's function, wavy line for 
the fluctuation propagator K, the shaded rectangle for the 
Cooperon C and shaded triangle for the vertex C/2ttut. 



the Fermi surface, and q is the momentum in the effec- 
tive dimension. Using Ea. (|lt)(l in a standard way, one 
obtains the "fluctuation propagator" i.e the impurity- 
averaged sum over the ladder diagrams corresponding to 
the electron-electron interaction in the Cooper channel, 



K^i^q) = In 



T 



■1> 



mi/2 



2ttT 

(11) 

where Q is the bosonic Matsubara frequency. The pole 
of Ea. Hll|) defines the boundary between the supercon- 
ducting and normal phases given by Eq. J2J) . At low tem- 
peratures, T <C a c o, the fluctuation propagator reduces 
to 

K-\n, q) = \n[{a q + |fi)|/2)/a c (T)]. (12) 

The fluctuation correction to the conductivity is ob- 
tained as usual from the Kubo formalism with the ap- 
propriate analytic continuation. The standard set of di- 
agrams constituting the AL, DOS and MT contributions 
is shown in Fig. 3. 

The AL contribution (Fig. 3a) can be expressed as a 
sum of two terms: 5<j al — Sa AL 
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where r 0q = (q/27rT)V>'[I /2 + {a q - zO/2)/2ttT]. The 
contribution 5<j al appears only at finite temperature 
while 5<72 L contains the contribution that survives even 
at T — 0. Note that Sa AL that describes the quantum 
regime, results from differentiating r q (— ifl) w.r.t f2 and 
is missed within the usual static approximation ( Q = 0) . 



The MT correction is given by 8a MT 
where 
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with uj + = uj + Q/2. For low temperatures, the contri- 
butions 6af IT (Fig. 3b) and 8a\ 1T (Fig. 3c) are of the 
same order and both need to be taken into account. On 
the contrary, at high temperatures, the diagram (Fig. 3c) 
having an extra Cooperon propagator is of a lower order. 

The DOS fluctuation correction is given by the expres- 
sion 6a DOS = 5a? os + 5a§ os where 
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The correction 8<j®° corresponds to diagrams shown in 
Fig. 3d, 3e while (5cr 2 DO ' S is given by diagrams from Fig. 
3f having an extra Cooperon propagator. 

Equipped with the analytical expressions for all the 
necessary diagrams, we now proceed to their evaluation 
beginning with the zero temperature case. One can see 
that apart from Sa AL , only the terms in the DOS and 
MT contributions that have a cth(f2) survive at T = 0. 
fntegrating over frequency f2 by parts, the T = answers 
for all contributions can be expressed in terms of just two 
functions defined by 



F 2 



2D 
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to yield 8a AL = 6F 2 , Sa uub = -(Fi + 3F 2 )/2, and 
Sa MT = ~3F 2 19\. Since the functions F 1 and F 2 are 
positive, we find that the AL correction is positive and 
the DOS correction is negative as expected while the MT 
correction which doesn't have a prescribed sign is nega- 
tive. Finally the fluctuation correction to the normal 
state conductivity given by Sa' ° tal =-(Fi~ 3F 2 )/2, af- 
ter integrating by parts results in Eq. J3J . 
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At nonzero temperatures, the terms having a factor 
l/sh(f2/2T) need to be included. One finds that the lead- 
ing contribution close to the QCP comes from 5af L and 
the evaluation of this term in different regimes leads to 
the results Eq.lQHJ) discussed earlier. 

Finally we would like to point out an interesting ex- 
perimental realization of QCP in a quantum wire, that 
of a hollow cylinder with thin walls |3J. In this case the 
pair-breaking parameter that can be easily obtained form 
the Usadel equation |2ij, reads 



a = D 



eH 

4c 



-An ■ 



— {rl + rl) 
c 



,ln(r 2 /ri) 



(21) 

where r% and ri are the inner and outer radii respectively 
and n is an arbitrary integer. For a thin cylinder (r\ w 
r<x « r) it reduces to 



a= (£>/2r 2 )($/$ - nf 



(22) 



where $o = nc/e and <f> is the flux enclosed by 
the cylinder, thereby rendering the classic Little-Park 
oscillations 21] of T c as can be seen from Eq.J2J). In- 
terestingly, for a cylinder with small enough radius, 
r < r c = yj Dj/AirT co , it is possible to push the T c down 
to zero at magnetic fields corresponding to half-integer 
fluxes <E> = $o(l/2 + n), as was experimentally observed 
in Ref. While the positive fluctuation contribution 
to conductivity that we would associate with the classi- 
cal regime was clearly observed in this experiment, the 
conductivity behavior expected for the quantum and in- 
termediate regimes were not reported so far. 

In conclusion, we have investigated the fluctuation cor- 
rection to the normal state conductivity in the vicinity 
of a parallel-field-induced QCP in dirty samples of re- 
duced dimensions, taking into account both, quantum 
and thermal fluctuations within the diagrammatic per- 
turbation theory. Our key finding is that there are 
three regimes that show a qualitatively different behavior 
ranging from quantum to classical. The particular tem- 
perature and field behavior of the conductivity is dic- 
tated by the choice of path in approaching the QCP 
while making the measurement. We have found that 
for a nanowire (or for a hollow cylinder) as well as for 
a thin film, the zero temperature conductivity correc- 
tion that also governs the quantum regime, is negative, 
which means that the quantum pairing-fluctuations in- 
crease the resistance to the charge flow. Our findings im- 
ply that experiments should detect a negative magnetore- 
sistance in the quantum regime. For detailed compari- 
son of our results for conductivity dependence on tem- 
perature and magnetic field with the experimental data, 
the weak localization!^ and Altshuler- Aronov^^ cor- 
rections must be subtracted from the experimental con- 
ductivity dependence. Inclusion of these corrections will 
not affect the predicted negative sign of the magnetore- 
sistance at low temperatures since the weak localization 



correction also results in a negative magnetoresistance 
while the Altshuler-Aronov correction does not depend 
on the magnetic field. 
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